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The  coherence  of  electron  spin  can  be  significantly  enhanced  by  locking  the  Overhauser  field  from 
nuclear  spins  using  the  nuclear  spin  preparation.  We  propose  a  theoretical  model  to  calculate  the 
long  time  dynamics  of  the  Overhauser  held  under  intrinsic  nuclear  spin  diffusion  in  a  quantum  dot. 
We  obtain  a  simplified  diffusion  equation  that  can  be  numerically  solved  and  show  quantitatively 
how  the  Knight  shift  and  the  electron-mediated  nuclear  spin  Hip-hop  affect  the  nuclear  spin  diffusion. 
The  results  explain  several  recent  experimental  observations,  where  the  decay  time  of  Overhauser 
Held  is  measured  under  different  conhgurations,  including  variation  of  the  external  magnetic  held, 
the  electron  spin  conhguration  in  a  double  dot,  and  the  initial  nuclear  spin  polarization  rate. 

PACS  numbers:  73.21.La,  76.20. +q,  76.60.-k,  85.35.Be 


I.  INTRODUCTION 

Electron  spins  in  single  quantum  dots  consist  of  one  of 
the  promising  systems  for  realization  of  quantum  com¬ 
putation  |Tj.  The  spin  state  of  a  single  electron  in  a 
quantum  dot  can  be  coherently  controlled  either  opti¬ 
cally  through  fast  laser  pulses  or  electronically  through 
tuning  of  gate  voltages  mini.  In  experiments, 
the  coherence  time  of  the  electron  spin  is  limited  by  its 
hyperfine  coupling  to  the  nuclear  spin  environment  in 
the  host  semiconductor  material.  The  coupling  causes 
spectral  diffusion  and  gives  a  typical  spin  decoherence 
time  T%  ~  15?rs  for  the  electron  spin  qubit  [3  E] .  This 
coherence  time  could  be  significantly  prolonged  with  ap¬ 
plication  of  spin  echo  or  other  dynamic  decoupling  tech¬ 
niques  [5].  However,  implementation  of  these  techniques 
requires  repeated  applications  of  many  laser  pulses.  Each 
pulse  inevitably  induces  some  noise  by  itself,  which  lim¬ 
its  the  practical  performance  of  the  dynamic  decoupling 
techniques  to  suppress  spin  noise  under  real  environ¬ 
ments. 

Another  technique  to  increase  the  coherence  time  for 
the  electron  spin  is  through  the  dynamic  nuclear-spin 
preparation  (DNP)  [9]  U0j  [TTJ  Q21  [131,  HU  [15],  which 
prepares  the  nuclear  spin  environment  into  certain  con¬ 
figurations.  The  nuclear  spins  in  this  configuration  col¬ 
lectively  generate  an  effective  magnetic  field  (the  Over¬ 
hauser  field)  on  the  electron  spin  with  small  fluctuation, 
significantly  reducing  the  spectral  diffusion  of  the  elec¬ 
tron  spin  qubit.  Although  nuclear  spins  can  be  polarized 
through  many  methods  such  as  optical  pumping  [5] ,  sub¬ 
stantial  reduction  in  the  fluctuation  of  the  Overhauser 
field  requires  almost  complete  polarization  of  the  nuclear 
spins  m,  which  is  hard  to  achieve  experimentally.  The 
recent  experiments,  however,  demonstrate  a  surprising 
feedback  mechanism  which  can  lock  the  Overhauser  field 
to  certain  values  without  significant  polarization  of  the 
nuclear  spin  environment  PU  HI  HU  HU-  The  Over¬ 
hauser  field  generated  from  this  locking  mechanism  has 
small  fluctuation,  which  effectively  increase  the  coherence 
time  T£  of  the  electron  spin  qubit  by  up  to  two  orders  of 


magnitudes.  The  locking  of  the  Overhauser  field  through 
the  DNP  process  has  been  seen  in  a  number  of  experi¬ 
mental  systems,  under  both  optical  and  electronic  gate 
control  [T21  IT3l  H4l  IT5]. 

To  enjoy  a  longer  spin  coherence  time,  it  is  desir¬ 
able  to  do  gate  experiments  on  the  electron  spins  under 
the  nuclear  spin  environment  prepared  with  the  DNP, 
which  produces  a  fixed  Overhauser  field.  An  important 
question  is  then  how  long  this  fixed  Overhauser  field 
can  survive  after  the  DNP  process  and  what  facts  de¬ 
termine/influence  the  relaxation  time  of  the  Overhauser 
field.  Recent  experiments  have  observed  that  the  pro¬ 
duced  Overhauser  field  can  endure  from  a  few  seconds  to 
a  few  minutes  and  the  relaxation  time  of  the  Overhauser 
field  depends  on  a  variety  of  experimental  parameters, 
such  as  the  applied  magnetic  field  EfflnanE],  the  elec¬ 
tron  spin  configuration  in  double  quantum  dots  m ,  and 
the  DNP  pump  time  [FT] . 

In  this  paper,  we  develop  a  quantitative  theory  to  cal¬ 
culate  the  relaxation  time  of  the  Overhauser  field  un¬ 
der  the  environment  of  quantum  dots,  and  provide  an 
explanation  to  the  experimental  observations  mentioned 
above  in  different  setups  and  with  control  of  different 
experimental  parameters.  The  relaxation  of  the  Over¬ 
hauser  field  is  caused  by  nuclear  spin  diffusion,  which 
has  been  well  studied  in  bulk  material  M  ED,  where 
diffusion  is  caused  simply  by  the  nuclear  dipole-dipole 
interaction.  In  a  quantum  dot,  however,  the  presence 
of  the  electron  spin  generates  several  new  effects.  First, 
the  electron  spin  can  mediate  the  diffusion  of  the  nuclear 
spins  through  a  virtual  hyperfine  process.  Second,  the  ef¬ 
fective  magnetic  field  from  the  electron  spin  produces  an 
inhomogeneous  Knight  shift  on  the  nuclear  spins  through 
the  hyperfine  coupling  which  suppresses  the  nuclear  spin 
diffusion.  The  change  of  the  nuclear  spin  diffusion  coef¬ 
ficient  by  the  Knight  shift  has  been  taken  into  account 
in  Ref.  [22],  but  without  consideration  of  the  electron 
mediated  nuclear  spin  diffusion.  A  recent  work  considers 
relaxation  of  the  Overhauser  field  due  to  the  electron  me¬ 
diated  nuclear  spin  diffusion,  but  without  consideration 
of  the  direct  nuclear  dipole-dipole  interaction  [23] .  In  this 
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case,  the  Overhauser  field  can  only  decay  by  less  than  1%, 
which  suggests  that  it  is  necessary  to  include  the  dipole 
interaction  in  the  long  time  dynamics.  A  quantitative 
theory  is  still  missing  to  our  knowledge  that  includes  a 
complete  description  of  all  the  competing  effects  men¬ 
tioned  above.  In  this  work,  we  take  into  account  all  these 
diffusion  mechanisms,  and  the  resulting  theory  provides 
an  explanation  of  the  recent  experimental  observations 
in  Refs.  mmmm- 

The  paper  is  arranged  as  follows:  in  Sec.  II  we  give 
a  formalism  to  describe  relaxation  of  the  Overhauser 
field  that  includes  contributions  from  the  nuclear  dipole- 
dipole  interaction,  the  electron  mediated  nuclear  spin  dif¬ 
fusion,  and  the  Knight  shift.  The  effective  nuclear  spin 
diffusion  equation  is  solved  numerically  to  determine  the 
relaxation  time  of  the  Overhauser  field.  In  Sec.  Ill,  we 
compare  the  theoretical  calculations  with  the  recent  ex¬ 
perimental  observations  under  control  of  different  experi¬ 
mental  parameters,  and  show  that  they  are  in  qualitative 
or  semi-quantitative  agreement.  We  summarize  our  re¬ 
sults  in  Sec.  IV  with  brief  discussions. 


II.  DECAY  OF  THE  OVERHAUSER  FIELD 
THROUGH  NUCLEAR  SPIN  DIFFUSION 


We  assume  that  an  external  magnetic  field  Bq  much 
larger  than  the  mean  value  and  variance  of  the  local 
Overhauser  field  generated  by  nuclei  is  applied  along  the 
z-direction  (perpendicular  to  the  quantum  dot  layer).  In 
this  case,  we  can  drop  the  nonsecular  terms  in  the  inter¬ 
action  Hamiltonian  |24|.  For  simplification,  we  consider 
only  one  species  of  nuclei  around  the  quantum  dot  elec¬ 
tron.  The  total  Hamiltonian  for  the  electron  and  nuclear 
spin  system,  including  both  the  Fermi  contact  hyperfine 
interaction  and  nuclear  dipole-dipole  interaction,  can  be 
written  as: 


H  —  He  +  Hn  +  Hen  +  Hnn,  (1) 

He  =  -gepBB0Sz,  (2) 

Hn  =  -gnmBoY.il  1  (3) 

i 

Hen  =  y A^ZIi  +  E  y(^+/r  +  s~i+),  (4) 

i  i 

Hnn  =  E  2iW>  ^  E  'V//,  >  (5) 
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Bij  =  ^(gnnN)2Rij3{l  ~  3cos2%),  (6) 

where  A;  denotes  the  hyperfine  coupling  between  the 
electron  and  nuclear  spin  at  site  i  with  spatial  coordi¬ 
nates  ( Xi,yi,Zi ).  Rij  is  the  distance  between  two  nuclei 
at  site  i,j.  9ij  is  the  angle  between  the  line  connecting 
sites  i,j  and  the  z  direction. 

We  note  that  for  B0  ranging  from  a  few  mT  to  a 
few  T,  the  electron  Zeeman  splitting  is  on  the  order  of 


1CT1  —  102GHz,  while  the  average  hyperfine  coupling  in 
most  quantum  dot  systems  is  on  the  order  of  MHz.  Thus 
we  can  adiabatically  eliminate  the  spin-flip  terms  in  the 
hyperfine  interaction  Hamiltonian  and  correspondingly 
modify  the  other  terms  in  the  Hamiltonian  as  [S3] : 
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where  we  have  introduced  an  electron-mediated  nuclear 
flip-flop  term  in  Hnn-  Since  we  are  interested  in  the  long 
time  dynamics  of  nuclear  spins,  we  can  completely  elimi¬ 
nate  the  electron  from  the  Hamiltonian  by  replacing  the 
constant  operator  S z  with  its  expectation  value.  We  find 
that  using  Sz  =  1/2  or  Sz  =  —1/2  will  yield  almost  the 
same  result  in  the  following  calculations.  Therefore  we 
can  set  Sz  =  1/2  for  simplicity  and  arrive  at  the  following 
effective  Hamiltonian  (neglecting  constant  terms): 
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Here,  the  term  proportional  to  A*  in  Hq  is  the  Knight 
shift  term.  For  this  Knight  shift,  we  have  neglected  the 
small  term  proportional  to  A2  in  Eq.  ^  as  it  is  domi¬ 
nated  by  the  Aj  term. 

The  expectation  value  for  z  component  of  the  nuclear 
spin  at  site  k  will  evolve  according  to  the  Schrodinger 
equation: 

^L=l-Tr{p{t)[H1Jzk}},  (13) 

where  p{t)  is  the  nuclear  spin  density  matrix  at  time  t , 
which  can  be  calculated  by  switching  to  the  interaction 
picture: 

p(t)=p(0)+l-j\p(t'),H1(t')}dt',  (14) 

with  H\[t)  =  exp(iHot/h)Hiexp(—iHot/h).  Further  cal¬ 
culation  yields  [2D]: 

^l=l-Tr{PmH1(t),Izk)} 

+  Q)  J*Tr{p(t-l')[H1,[H1(l'))lZ\]}dt?  (15) 

We  assume  the  nuclear  spin  (with  spin-/)  density  ma¬ 
trix  in  a  product  state  of  the  following  form: 

M  =  nit)  =  ^  (16) 
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Such  an  approximation  is  valid  for  the  initial  state  of  nu¬ 
clear  spin  when  correlations  and  transverse  coherence  are 
negligible  in  its  equilibrium  configuration.  Off-diagonal 
terms  may  appear  in  p(t)  for  t  >  0,  but  it  can  be  shown 
that  they  only  have  a  minor  contribution  to  the  evolution 
of  (I|(t))  compared  to  the  diagonal  part  expressed  in  Eq. 
( 16 1  (see  Ref.  [20]). 
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By  using  the  explicit  form  of  the  Hamiltonian  [Eq.  ( 10 
and  density  matrix  [Eq.  ©]>  we  can  reduce  Eq. 


(151  to: 


=  XVifc«I?(i)>-<Jfc(t)>)  (17) 

i^k 

Wkl  =  Tr{(P)2}  l  [Hi(t  ~ 


where  Wki  has  a  clear  physical  meaning  as  the  flip-flop 
rate  between  nuclear  spins  at  site  i  and  k. 

For  a  2D  InAs/GaAs  quantum  dot,  we  take  As  nuclei 
(J  =  3/2)  as  an  example  for  further  calculation.  The 
parameter  Wki  can  be  analytically  calculated  when  ap¬ 
proximating  the  integration  upper  limit  in  the  above  ex¬ 
pression  for  Wki  from  t  to  infinity  [22], 
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The  hyperfine  coupling  rate  A,  is  proportional  to  square 
of  the  electron  wave  function  in  a  quantum  dot.  And  in 
the  following  calculation  we  assume  the  dot  potential  is 
like  a  square  well  in  the  ^-direction  and  the  electron  wave 
function  takes  a  Gaussian  shape  in  the  cc,  y-plane.  The 
hyperfine  coupling  rate  Ai  can  then  be  written  as  A,  = 
A0 cos2 (ir Zi / Zo)  exp[- (x? +y?)/l%],  where  (xi:yi,Zi)  are 
spatial  coordinate  of  the  site  i.  Iq  and  zq  are,  respectively, 
the  Fock-Darwin  radius  and  thickness  of  the  quantum 
dot,  and  Aq  is  the  hyperfine  coupling  for  the  nuclear  spin 
at  the  origin  where  the  electron  locates. 

We  assume  the  nuclei  spins  follow  a  diffusion  process, 
where  the  flip-flop  rate  for  two  distant  sites  is  negligible. 
This  approximation  is  justified  by  the  fact  that  the  coeffi¬ 
cient  generally  decays  to  zero  fast  as  Rik  increases.  If 
we  treat  (Ik{t))  as  a  continuous  function  of  spatial  vari¬ 
able  xa  (a  =  x,y,z).  We  can  then  carry  out  a  Taylor 


expansion  of  ( Iz{t ))  for  site  i  around  site  k: 
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1  d  (Ik  (t))  _  p  _  p  ,  _  C22) 

+  2dxadxf>{k  iKk  {  ) 

where  Einstein’s  summation  convention  is  implied  for 
spatial  index  a,  (3.  Substituting  this  into  Eq.  (|T7|)  and 
noting  that  the  summation  of  the  first  order  derivative 
term  over  all  sites  vanishes  due  to  the  lattice  symmetry 
[25].  we  have: 
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The  ~k  notation  above  means  summarization  over  the 


k  - 


sites  near  k.  Define  coefficient  Da/3  =  Y2i*ik  W/ 
x?)(xk  —  xi)/2  and  similarly  note  that  for  a  /  /3  the 
summation  over  all  sites  vanishes,  we  have: 


mi 

dt 


cB 

=  (Dxx-—  +  Dyv-— 
1  dx 2  dy 2 


Dzz^){m)  (24) 


Eq.  (24 1  is  a  3D  anisotropic  diffusion  equation  with 


spatially  varying  diffusion  coefficients  (as  Aikl  Bik,  Wik 
all  depend  on  the  spatial  coordinates),  which  is  not  easy 
to  solve.  To  further  simplify  it,  we  note  that  to  obtain 
the  major  feature  for  the  full  time  dynamics  of  the  Over- 
hauser  field  ( hz(t ))  =  J2k  Ak{Ik{t)),  it  is  reasonable  to 
first  ignore  the  diffusion  in  the  z  direction  as  the  quan¬ 
tum  dot  layer  is  usually  a  few  nm  thick  and  chemical 
or  structural  mismatch  in  adjacent  layers  may  strongly 
suppress  diffusion  in  the  z  direction  |I8I.  In  addition, 
from  symmetry  in  the  2D  x,  y-plane,  we  expect  to  have 
Dxx  Ki  Dvv  and  thus  define  an  average  2D  diffusion  co¬ 
efficient  D(x,y )  =  Wik[(xk  -  Xi )2  +  (yk  -  2/i)2]/4. 
Now  we  have  a  simplified  2D  diffusion  equation: 


d{I} 


dT  =  D^dx* 


d2  d2 


dy 


(25) 


The  above  partial  differential  equation  can  be  effective 
solved  using  the  finite  element  method  by  coarse  graining 
a  large  number  of  nuclear  spin  sites  to  a  small  number 
of  mesh  nodes.  But  before  solving  Eq.  (25 1,  we  would 


like  to  have  some  discussion  about  the  diffusion  coeffi¬ 
cient  D{x,y).  For  i,|/>  l0,  the  role  of  electron  can  be 
neglected  and  a  numerical  calculation  of  the  above  diffu¬ 
sion  coefficient  gives  a  uniform  background  value  D  «  7 
nm2 /s,  which  is  consistent  with  the  previous  theoretical 
and  experimental  study  of  the  diffusion  coefficient  in  the 
bulk  material  I2H1I2HE2]-  In  our  calculation,  we  set  the 
quantum  dot  parameters  as  Iq  =  30  nm,  Zq  =  10  nm, 
A0  =  1  yeV  ~  1.5  MHz,  J2k^k  ~  80  yeV,  the  lat¬ 
tice  constant  ao  =  0.563  nm,  and  the  number  of  nuclei 
TV  a:  9  x  105,  according  to  the  typical  experimental  values 

Usmans]. 

Within  the  range  of  Fock-Darwin  radius  Iq,  however, 
the  presence  of  the  quantum  dot  electron  will  change  the 
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diffusion  coefficient  through  two  competing  mechanisms: 
on  the  one  hand,  the  presence  of  the  confined  electron 
generates  an  inhomogeneous  Knight  shift  [B|,  which  lifts 
the  degeneracy  of  the  nuclear  Zeeman  energy  for  differ¬ 
ent  nuclei  and  prevents  the  spin  flip-flop;  on  the  other 
hand,  electron  mediated  nuclear  spin  flip-flop  enhances 
the  nuclear  spin  diffusion  from  the  center  to  the  edge  of 
the  dot. 


Our  numerical  simulation  shows  that  whether  one 
mechanism  dominates  the  other  is  largely  determined  by 
the  external  magnetic  field  Bq.  Fig.  m  show  the  diffu¬ 
sion  coefficient  D(x,y)  under  Bq  =  0.2 T  and  Bq  =  2 T. 
We  can  see  that  the  electron  mediated  flip-flop  greatly 
enhances  the  nuclear  spin  diffusion  near  the  center  of  the 
dot  under  a  small  magnetic  field  (with  a  sharp  peak  in 
D(x,y)),  while  for  a  large  magnetic  field,  this  effect  is 
negligible  compared  to  the  Knight  shift  which  suppresses 
the  nuclear  spin  diffusion  (with  a  wide  dip  in  D(x,y)). 
The  difference  can  be  easily  explained  from  the  effective 
Hamiltonian  [Eq.  ( To]|l2|]:  the  electron  mediated  flip- 
flop  term  is  inversely  proportional  to  Bq  while  the  Knight 
shift  term  is  independent  of  B0.  We  note  that  the  rea¬ 
son  why  we  have  a  narrower  peak  than  the  dip  is  due  to 
the  fact  that  the  Knight  shift  term  is  proportional  to  the 
hyperfine  coupling  (with  a  Gaussian  distribution)  while 
the  Electron  mediated  flip-flop  term  is  proportional  to 
the  product  of  two  nuclei’s  hyperfine  coupling  rates.  We 
also  note  that  the  2D  diffusion  coefficient  D(x,y)  in  Fig. 
[l]|2]  does  not  have  azimuthal  symmetry  since  in  our  cal¬ 
culation  we  assume  the  nuclear  spins  in  a  square  lattice 
which  has  no  azimuthal  symmetry. 


III.  COMPARISON  WITH  EXPERIMENTS 


To  compare  with  experiments,  we  numerically  solve 
the  diffusion  equation  [Eq.  (25 1]  under  certain  initial 
and  boundary  conditions.  For  the  initial  condition,  the 
nuclear  spins  are  partially  polarized  through  the  DNP 
process  from  the  hyperfine  interaction  with  the  elec¬ 
tron  spin  [H  EH  M  CD2-  It  is  reasonable  to  expect 
that  right  after  the  DNP  process,  the  polarization  dis¬ 
tribution  (/|)  is  proportional  to  the  hyperfine  interac¬ 
tion  rate.  So  in  the  following  calculation,  we  assume 
(ID  cx  Ak  oc  exp[—  (x'j  +  iiD/Iq]  at  t  =  0  for  solving 
the  diffusion  equation  [Eq.  (j25j[.  We  can  take  the  natu¬ 
ral  boundary  condition  where  (Ik)  approaches  zero  when 
the  radius  goes  to  infinity.  However,  in  numerical  calcu¬ 
lation,  we  have  to  take  a  finite  radius.  To  make  the  spin 
diffusion  possible,  this  finite  radius  has  to  be  significantly 
larger  than  the  radius  of  the  size  Iq  of  the  initial  electron 
wave  packet.  In  the  calculation,  we  typically  take  the 
integration  size  about  300  nm  (10  times  of  Iq)  so  that 
the  number  of  total  nuclear  spins  in  this  area  is  about 
100  times  the  number  of  initially  partially  polarized  nu¬ 
clear  spins  covered  by  the  electron  wave  packet.  With 
this  choice,  we  expect  the  cutoff  error  should  be  small 
(at  a  percent  level) . 
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FIG.  1:  Diffusion  Coefficient  D(x,y)  under  Bq  =  0.2 T.  The 
narrow  high  peak  at  the  center  of  the  dot  is  due  to  electron- 
mediated  nuclear  spin  flip-flop,  and  the  wide  dip  is  due  to 
inhomogeneous  Knight  shift. 


FIG.  2:  Diffusion  Coefficient  D(x,y)  under  Bo  =  2 T.  The  in¬ 
homogeneous  Knight  shift  dominates  in  this  case,  so  diffusion 
is  generally  suppressed  within  the  Fock-Darwin  radius. 


First,  to  compare  with  the  experiments  in  Ref.  m  ,  we 
calculate  the  relaxation  of  the  Overhauser  field  hz(t)  = 
J2k  Ak(Ik(t)}  under  different  electron  states  and  different 
values  of  the  external  magnetic  field  Bq,  and  the  result 
is  shown  in  Fig.  [3j  We  note  that  for  the  double  quan¬ 
tum  dot  system  in  Ref.  m  if  the  electron  stays  in  the 
(2,0)  singlet  state,  the  electron  spin  has  Sz  =  0,  with 
basically  no  influence  on  the  nuclear  spin  diffusion.  In 
this  case,  the  nuclear  spin  diffusion  is  governed  by  the 
intrinsic  nuclear  dipole-dipole  interaction.  However,  for 
the  electron  in  the  (1,1)  state,  with  the  magnetic  field 
in  the  range  of  tens  of  mT  as  in  this  experiment,  the 
electron  mediated  spin  diffusion  dominates  the  Knight 
shift  and  it  accelerates  the  nuclear  spin  relaxation.  One 
can  see  two  effects  from  Fig.  [3]  with  the  magnetic  field 
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Time  (s) 


FIG.  3:  Decay  of  the  Overhauser  field  under  various  small 
magnetic  fields.  The  solid  line  refers  to  the  case  with  electron 
staying  in  the  (2,0)  singlet  state  where  electron  plays  no  role 
in  nuclear  spin  diffusion. 


FIG.  5:  Overhauser  Decay  of  the  Overhauser  field  under  vari¬ 
ous  initial  distributions  of  the  nuclear  polarization.  The  solid 
line  refers  to  the  case  with  a  long  DNP  pump  time  that  gives 
a  Gaussian  distribution  with  the  size  characterized  by  the 
Fock-Darwin  radius  lo .  Other  lines  correspond  to  narrower 
polarization  distribution  characterized  by  a  Gaussian  with  its 
size  ?'o  <  lo- 


FIG.  4:  Decay  of  the  Overhauser  field  under  a  strong  magnetic 
field.  The  solid  line  shows  the  dynamics  of  the  Overhauser 
field  under  B  =  2 T  and  the  dotted  line  is  a  fit  by  using  a 
constant  diffusion  coefficient. 


Bq  =  10  mT  or  20  mT:  (1)  electrons  in  the  (1,1)  state  will 
speed  up  the  decay  of  the  Overhauser  field  compared  to 
electrons  in  the  (2,0)  state;  (2)  a  smaller  magnetic  field 
gives  a  faster  decay  of  the  Overhauser  field.  Both  of  these 
effects  agree  with  the  experimental  observations  in  Ref. 
sa.  and  the  decay  time  scale  is  also  consistent  with  what 
one  measures  from  the  experiments  in  term  of  the  order 
of  magnitude. 

With  a  much  larger  magnetic  field  (say,  Bq  =  2 T,  as  in 
experiments  in  Ref.  |18j).  the  electron  mediated  nuclear 
spin  diffusion  is  suppressed,  and  the  Knight  shift  plays  a 
more  important  role.  The  Knight  shift  suppresses  the  nu¬ 
clear  spin  diffusion,  and  can  prolong  the  relaxation  time 


of  the  Overhauser  field  to  make  it  significantly  larger  than 
the  relaxation  time  in  the  bulk  material.  Fig.  [4]shows  de¬ 
cay  of  the  Overhauser  field  in  this  case,  and  we  can  fit  the 
curve  with  an  effective  constant  diffusion  coefficient  with 
its  value  about  Deff  «  0.7 nm2/s.  Compared  with  the 
diffusion  coefficient  in  the  bulk  material  ( D  ss  7 nm2 /s), 
the  effective  diffusion  here  is  suppressed  by  a  factor  of  10 
under  a  strong  magnetic  field  B0.  Experiments  done  in 
Ref.  [T5]  yield  an  effective  diffusion  coefficient  50  times 
smaller  than  the  value  in  the  bulk  system.  While  the 
suppression  there  could  have  contribution  from  other  fac¬ 
tors,  such  as  inhomogeneity  of  lattice  constant  there,  we 
believe  that  the  Knight-shift  induced  suppression  plays 
a  large  role  in  this  experiment  under  a  strong  magnetic 
field. 

The  experiment  in  Ref.  [T6]  studies  relaxation  of  the 
Overhauser  field  under  different  pumping  time  for  the 
DNP  process.  With  a  shorter  DNP  pumping  time,  the 
nuclear  spin  polarization  may  have  a  narrower  distribu¬ 
tion  in  space  [IB].  Although  we  do  not  know  the  exact 
distribution  of  the  nuclear  spin  polarization  from  a  short 
DNP  pump,  to  model  this  effect  qualitatively,  we  simply 
assume  that  this  distribution  (/|)  is  still  a  Gaussian  but 
with  its  radius  ro  <  lo-  Taking  this  (/f )  as  the  initial  con¬ 
dition,  we  can  calculate  relaxation  of  the  corresponding 
Overhauser  field  from  the  diffusion  equation  [Eq.  ( 25 1] , 
and  the  result  is  shown  in  Fig.  [5]  The  result  indicates 
that  a  narrower  distribution  of  initial  nuclear  spin  polar¬ 
ization  leads  to  a  faster  decay  of  the  Overhauser  field, 
which  is  consistent  with  the  experimental  result  in  Ref. 
[16] .  This  effect  can  be  explained  by  noting  that  the  dif¬ 
fusion  is  much  stronger  near  the  center  of  the  dot  due  to 
the  electron  mediated  diffusion  peak  (see  Fig.  [TJ ,  which 
results  in  a  shorter  relaxation  time  of  the  Overhauser 
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field  if  the  initial  polarization  is  more  concentrated  near 
the  dot  center. 


IV.  SUMMARY  AND  DISCUSSION 

In  summary,  we  have  established  an  effective  method 
for  calculating  the  long  time  dynamics  of  the  Overhauser 
field  under  nuclear  spin  diffusion  and  shown  that  the  con¬ 
fined  electron  in  a  quantum  dot  can  both  enhance  decay 
of  the  Overhauser  field  by  mediating  nuclear  spin  flip-flop 
and  suppress  nuclear  spin  diffusion  via  inhomogeneous 
Knight  shift.  Which  effect  dominates  depends  critically 
on  the  magnitude  of  the  external  magnetic  field.  We  sim¬ 
ulate  the  relaxation  of  Overhauser  field  under  different 
electron  spin  configuration,  external  magnetic  field,  and 
initial  nuclear  polarization  distribution,  and  the  results 
agree  qualitatively  with  a  series  of  recent  experimental 
observations. 

For  the  purpose  of  maintaining  the  Overhauser  field 
generated  by  the  DNP  process  as  long  as  possible,  ap¬ 
plying  a  large  magnetic  field  turns  out  to  be  the  most 
effective  way  since  a  large  magnetic  field  suppresses  the 
electron  mediated  nuclear  spin  diffusion.  The  intrinsic 
nuclear  spin  diffusion  from  the  nuclear  dipole-dipole  in¬ 
teraction  is  also  suppressed  by  the  inhomogeneous  Knight 
shift.  These  two  kinds  of  suppression,  combined  together, 


leads  to  a  long  relaxation  time  of  the  Overhauser  field. 

For  calculations  in  this  paper,  we  focus  on  the  relax¬ 
ation  dynamics  of  the  expectation  value  of  the  Over¬ 
hauser  field,  since  this  is  the  quantity  that  has  been 
measured  in  several  recent  experiments.  Similar  meth¬ 
ods  could  also  apply  to  calculation  of  relaxation  of  the 
variance  of  the  Overhauser  field,  and  in  term  of  the  time 
scale,  the  relaxation  time  for  variance  is  basically  the 
same  as  the  relaxation  time  for  the  expectation  value. 

During  the  DNP  process,  nuclear  spin  diffusion  also 
takes  place.  The  final  distribution  of  the  nuclear  spin  po¬ 
larization  and  the  limit  on  its  fluctuation  may  depend  on 
a  balance  between  the  DNP  pumping  cycle  and  the  con¬ 
tinuous  nuclear  spin  diffusion  process  m  ei  mi  asi  Ei- 
To  understand  this  balance,  we  need  a  more  detailed  un¬ 
derstanding  of  the  dynamics  for  the  DNP  process  and  its 
dependence  on  various  experimental  control  parameters. 
This  is  an  interesting  topic  for  further  investigation. 


Acknowledgments 

We  thank  Duncan  Steel  for  helpful  discussions.  This 
work  is  supported  by  the  ARO  MURI  program,  the 
IARPA  grants,  the  DARPA  OLE  Program,  and  the 
AFOSR  MURI  program. 


[1]  D.  Loss  and  D.  P.  DiVincenzo,  Phys.  Rev.  A,  57,  120 
(1998). 

[2]  D.  D.  Awschalom,  et  al,  Semiconductor  Spintronics  and 
Quantum  Computation  (Springer,  Berlin,  2002). 

[3]  J.  M.  Taylor,  et  al,  Nature  Physics,  1,  177  (2005). 

[4]  J.  R.  Petta,  et  al,  Science,  309,  2180  (2005). 

[5]  X.-D.  Xu,  et  al,  Phys.  Rev.  Lett.,  99,  097401  (2007). 

[6]  A.  V.  Khaetskii,  D.  Loss,  and  L.  Glazman,  Phys.  Rev. 
Lett.  88,  186802  (2002). 

[7]  R.  Hanson,  et  al,  Rev.  Mod.  Phys.,  79,  1217  (2007). 

[8]  L.  Viola,  E.  Knill,  and  S.  Lloyd,  Phys.  Rev.  Lett.  82,  2417 
(1999);  L.-M.  Duan  and  G.-C.  Guo,  quant-ph/9807072, 
Phys.  Lett.  A  261,  139  (1999);  P.  Zanardi,  Phys.  Lett. 
A  258,  77  (1999);  H.  K.  Ng,  D.  A.  Lidar  and  J.  Preskill, 
arXiv:0911.3202  and  refs,  therein. 

[9]  A.  Imamoglu,  E.  Knill,  L.  Tian,  P.  Zoller,  Phys.  Rev. 
Lett.,  91,  017402  (2003). 

[10]  W.  M.  Witzel  and  S.  Das  Sarma,  Phys.  Rev.  B,  77, 
165319  (2008). 

[11]  J.  R.  Petta,  et  al,  Phys.  Rev.  Lett.  100,  067601  (2008). 

[12]  D.  J.  Reilly,  et  al.  Science,  321,  817  (2008). 

[13]  Xiaodong  Xu,  et  al,  Nature,  459,  1105  (2009). 


[14]  C.  Latta,  et  al,  Nature  Physics,  5,  758  (2009). 

[15]  I.  Vink,  et  al,  Nature  Physics,  5,  764  (2009). 

[16]  P.  Maletinsky,  A.  Badolato,  and  A.  Imamoglu,  Phys. 
Rev.  Lett.,  99,  056804,  (2007). 

[17]  D.  J.  Reilly,  et  al,  arXiv:0803.3082 

[18]  A.  E.  Nikolaenko,  et  al,  Phys.  Rev.  B.,  79,  081303(R), 
(2009). 

[19]  M.  N.  Makhonin,  et  al,  Phys.  Rev.  B  77,  125307  (2008) 

[20]  I.  J.  Lowe  and  S.  Gade,  Phys.  Rev.,  156,  817  (1967). 

[21]  D.  Pagnet,  Phys.  Rev.  B.,  25,  4444  (1982). 

[22]  C.  Deng  and  X.  Hu,  Phys.  Rev.  B.,  72,  165333  (2005). 

[23]  D.  Klauser,  W.  A.  Coish,  and  D.  Loss,  Phys.  Rev.  B., 
78,  205301  (2008). 

[24]  C.  P.  Slichter,  Principles  of  Magnetic  Resonance  (Harper 
and  Row  Inc.,  New  York,  1963). 

[25]  The  hyperfine  coupling  rate  Ai,  although  inhomogeneous 
over  a  larger  scale,  only  changes  slightly  from  one  nuclear 
spin  to  its  nearby  spins,  so  in  summation  over  the  nearby 
lattice  sites,  we  still  have  an  approximate  lattice  symme¬ 
try  for  the  jump  rate  Wik- 


